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Abstract 

We present the results of afterpulse measurements performed as qualification test for 473 
inner detector photomultipliers of the Double Chooz experiment. The measurements include 
the determination of a total afterpulse occurrence probability as well as an average time distri- 
bution of these pulses. Additionally, more detailed measurements with different light sources 
and simultaneous charge and timing measurements were performed with a few photomultipli- 
ers to allow a more detailed understanding of the effect. The results of all measurements are 
presented and discussed. 

1 Introduction 

The reactor antineutrino disappearance experiment Double Chooz is designed to determine the 
leptonic mixing angle #13 [1J. A disappearance of reactor v e with a best fit value of 

sin 2 26»i 3 = 0.109 ± 0.030 (stat) ± 0.025 (sys) (1) 

was already presented in (2j [3] . 

The experimental design consists of two identical detectors with a 10.3 m 3 target each. The target 
material is a gadolinium doped liquid scintillator in which neutrinos are detected by inverse beta 
decay. To detect the scintillator light each target is surrounded by 390 photomultiplier tubes 
(PMT), which are 10 inch, hemispherical PMTs from Hamamatsu (type R7081MOD-ASSY) with 
a bialcali photokathode (Sb-K-Cs). For the two detectors a total of 780 PMTs + 23 spares were 
provided, about one half by German and the other half by Japanese groups. 

For a preselection and validation of the specifications of the 803 target PMTs two test setups were 
built: (A) at the "Max-Planck-Institut fur Kernphysik" (MPIK) in Heidelberg in collaboration 
with the "RWTH Aachen University" @] and (B) at the "Tokyo Institute of Technology" [5]. This 
article focuses on afterpulse measurements that have been performed with setup A. 



* Corresponding authors. 
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For 473 PMTs the probability of the occurrence of afterpulses has been measured. Furthermore 
the average arrival time distribution of the afterpulses and their amplitudes were determined. Two 
different analysis methods of the same data set will be presented in this article. Additionally more 
detailed investigations have been done including a combination of correlated charge and timing 
measurements of the afterpulses. Different light sources and hardware components were used to 
validate the consistency of the different measurements. The more detailed measurements were 
done only for a few selected PMTs. 

2 Afterpulse probability 

Photoelectrons created at the photocathode of a photomultiplier (PMT) may cause ionization of 
the residual gas between the photocathode and the first dynode. Due to the electric field the 
positive ions travel back to the photocathode and create secondary electrons, which result in after- 
pulses. The photoelectrons of the afterpulses may themselves create afterpulses, which are called 
afterpulses of higher order. 

In this article the expected number of afterpulses £ induced by a single photoelectron (SPE) is 
determined, including also higher order afterpulses. Furthermore, the afterpulse probability "P , 
considering only first order afterpulses, is calculated. We define "P to be the probability of mea- 
suring at least one afterpulse induced by an initial SPE. The relation between £ and V will be 
discussed at the end of this section. 

Multiple ionizations induced by one single photoelectron are treated as independent processes. 
This leads to the conclusion that the probability P^(n,k) of creating k first order afterpulses 
induced by n initial photoelectrons obeys Poissonian statistics: 

P>,fc) = fe^e-«-<*. (2) 

Given that each ionization process results in an afterpulse, fi represents the average number of 
ionization processes per photoelectron. Based on this equation the probability P ap (n) of measuring 
at least one afterpulse for n initial photoelectrons is given by 

P ap (n)=P ll (n,k> 1) = 1 - e-"-" . (3) 

We define V :— _P ap (l), so equation [3] can be expressed as 

P ap (n) = l-(1-P)». (4) 

To calculate the expected number of afterpulses we start with the expected value of first order 
afterpulses _Ei st (n) induced by n initial photoelectrons. Using equation [2] we find: 

oo 

Ei st {n) = • k = n ■ /x . (5) 

fe=0 

In order to include second order afterpulses it is necessary to consider that afterpulses can contain 
more than one photoelectron. According to our model, the afterpulse charge £ is thought to be 
independent of the number of initial photoelectrons n and the afterpulse order. We therefore 
introduce the expected number of photoelectrons contained by an afterpulse £ = X)fco(^ ' 6e)> 
while g e represents the time-averaged afterpulse charge density histogram. The expected number 
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of second order afterpulses can now be written as 

OO OO CO 

k=a k>=oe=o 

oo 

= nfx ■ £fi ■ Q ( = n/i 2 - J. (6) 

For the number of expected afterpulses -B,(n) of the i-th afterpulse order we can then conclude: 

Ei(n) = Ei_ x {n) ■ n ■ 1 = nfi ■ (jjl ■ If' 1 . (7) 
Summation of all orders leads to a total number of expected afterpulses of 

OO 

E(n) = n-»-Y,(»-ey- 1 = 1 ^ Ji 1 = n-E, (8) 
i=i I H't 

while S := E(\) is defined as expected number of afterpulses induced by a single initial photoelec- 
tron. Combining equation [3] and l8].E(n) can be expressed as a function of P&p(ri): 

E{n) = - "-1°(1-^p("))_ . (9) 
n + ln(l - P ap (n)) • I 

In our measurement we will record the occurrence of afterpulses in the time range up to 16 |J.s, as 
the selection specification of the afterpulse probability was defined for this interval. In section [3. 2. 3 1 
we will show that first order afterpulses occur only in this time range and thus the measurement of 
fi will not be affected by the measurement's constraint. However, afterpulses of higher order may 
occur later than 16 \x& after the initial pulse and will not be detected in the measurement. In order 
to account for this effect, we separate equation [8] in first order and higher order afterpulses E^: 

s = — L = /x + £hj (io) 
l — fj, ■ i 

E h = ==H-i-S. (11) 

1 — (J, ■ I 

We define £ A and E^ to be the expected number of afterpulses within the time range from to 
16\is. Since only a fraction of higher order afterpulses 

/ = 9s (12) 

will be detected, we redefine equation [TUj 

S-=, + E h .f= ^±^l . (13) 
1 — fj, ■ £ 

In our analyses we will determine the discrete temporal distributions \ii and Si for the expected 
number of first and first plus higher order afterpulses in a time bin i. Based on equation |11| the 
temporal distribution of the expected number of higher order afterpulses E^ i can be calculated 
via 

i 

E hti = 1 -^Ek ■ m-k - (14) 

fe=0 
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Figure 1: Scheme of the experimental setup. 



The fraction / is then computed via the expression 



E h 



E h . 



(15) 



which is independent of I. 

3 Timing measurements with 473 PMTs 

Timing measurements of 473 PMTs were performed to evaluate the time-dependent occurrence of 
afterpulses. Additionally, the total afterpulse occurrence for the time range from to 16 |J.s after an 
initial signal was determined to validate the PMT specification for the Double Chooz experiment. 
The selection specification demaned a total afterpulse probability of less than 10 %. 
In this section we will present two methods to analyse the same afterpulse data set. At first 
the expected value of afterpulses £ produced per SPE event is calculated. The second analysis 
method will bring forth the expected value of first order afterpulses /i, from which the afterpulse 
probability V can be derived. Moreover, the average number of photoelectrons £ per afterpulse 



can be determined using £ and "P (see section 3.2.3). 
3.1 Experimental setup 

The experimental setup of the timing measurements is sketched in figure ^ It is a fraction of the 
complete experimental setup which was build at MPIK to allow for the qualification tests of 473 
PMTs. For more detailed information about the complete setup and all performed measurements 
the reader may refer to [3]. 

The PMTs were mounted on a rack system inside a light tight and clcctromagnetically shielded 
Faraday room. It allowed the calibration of 30 PMTs in parallel. The data acquisition system as 
well as the power supply was installed outside the Faraday room. A splitter box (a module designed 
and built for the DC experiment by the CIEMAT group, Madrid) was responsible for decoupling 
the signal from the supply voltage. Each box was connected to one channel of a high voltage power 
supply (SY 2527 Universal Multichannel Power Supply System, CAEN) using a positive polarity 
to operate the PMTs. The signal outputs were connected via a lOx amplifier (sixteen channel 
amplifier, Model 776, Philips Scientific) to the measurement devices. 



For the afterpulse measurement two trigger boards (TB), which were developed for the Dou- 
ble Chooz trigger system [BJ, have been used. These boards contain 18 analogue input channels 
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with two discriminators each, a sum signal line of all input channels with four discriminators, a 
rate counter for each discriminator line, a FIFO memory, which is capable of storing 126 events 
simultaneously, and a logical processing unit for the discriminator signals. The information which 
is stored in the FIFO for each event contains the rate counters, the discriminator's status of the 
triggered event as well as the time difference to the previous event with a resolution of 16 ns. The 
logical processing unit is freely programmable and based on logical OR and AND operations. For 
the measurement only one discriminator for each analogue input channel was used. Each channel 
whose discriminator produces a signal creates an internal event. Storing the discriminator status 
makes it possible to determine the corresponding channel. For the calibration of 30 PMTs two TBs 
where used, while each was connected to 15 PMTs. The discriminator for each PMT channel was 
set to 25 % of a single photoelectron pulse (SPE pulse) . The TBs are able to generate a NIM-based 
logic signal, triggered by an external or internal (via software) signal. This signal is used to trigger 
the light source and its arrival time was measured as strobe in a free channel of each TB. 

As light source the beam of a 438 nm picosecond injection laser (PiLas) from Advanced Laser 
Diode Systems was sent through a diffuser to a perforated plate with 30 optical fibers connected. 
The other end of the optical fibers were positioned in front of the PMT, illuminating the full cath- 
ode. In order to yield SPE events as signal the laser intensity was adjusted to a level where only 
about one out of 10 trigger signals led to a PMT signal. The average signal-to-trigger ratio was 
Rt = 9.95 % ± 2.98 %, while the error was derived from the variation in the ratios of all PMTs. It 
was not possible to adjust the light intensity at each PMT individually, as 30 PMTs were illumi- 
nated by the same laser beam. 

The applied high voltages had been calibrated previously, so a gain of 10 7 was adjusted at each 
PMTg|. 

During each measurement one TB released a signal to trigger the laser. The strobe signal to 
each TB defines the start time t for one measurement cycle. For all subsequent signals passing 
the discriminator threshold the time difference to the preceding signal was stored up to a maximum 
time difference of 20 us. The average time difference between the start time to and the detection 
of the laser signal is 0.7 \is. Hence, the collection of all subsequent signals up to 16 |os after the 
laser signal is ensured. For one measurement typically 500 000 trigger signals were sent, which 
corresponded to an average number of 50 000 SPE signals for each PMT. 

3.2 Analyses and results 

The stored TB data contained all arrival time differences between all incoming signals. As a first 
step the time differences between the trigger signal t and all incoming PMT signals are plotted 
in figure [2j Based on the TB resolution a bin width of At = 16 ns was used for all following 
distributions. 

The raw distribution is dominated by a prominent laser signal peak. The peak maximum defines 
the mean arrival time of the laser signal. The intrinsic arrival time variation is dominated by the 
PMTs transit time spread of typically 2.8 ns [7] . Due to the TB resolution the measured arrival time 
distribution is broadened, which leads to a FWHM of the laser signal peak of 16 -32 ns. Together 
with a contamination of pre- and late-pulses, arriving in a time range between —30 ns and 80 ns 
around the main signal, the total width of the laser signal peak is about 150 ns. The first signal 
within this time window is classified as the initial laser signal. If an initial signal is detected, all 
subsequent signals up to 16 \is after this initial signal are interpret as possible afterpulses and stored 
in a histogram Si . All signals without an initial signal are treated as a distribution containing only 
dark noise signals and yields in a binned form the background histogram Bi. 
With a signal-to-trigger ratio i?x of about 10 %, only one out of 10 measurement cycles contains an 
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Figure 2: (a) Distribution of measured time differences between detected PMT signals and the strobe signal, (b) 
Zoom of figure (a). 



initial signal plus afterpulses, whereas the remaining 90% contain only dark noise signals. Hence, 
the dark noise contribution has to be rescaled by i?x to be suitable for a dark noise correction in 
the initial signal distribution. The number of cycles with and without an initial signal determine 
Rt, which is calculated for each PMT individually. 

The goal of this analysis is to evaluate afterpulse effects for SPE events as the afterpulse probability 
increases with larger number of initial photoelectrons (as shown in equation Bj. However, the data 
sample is slightly contaminated by multi-photoelectron events (NPE events). In order to correct 
this effect by means of equation [4] and [HJ the average number of created photoelectrons n per initial 
signal has to be computed: 

The probability p pe (n > 1) of creating at least one photoelectron for a released trigger signal is 
Poisson distributed, characterized by the average number of created photoelectrons per trigger 
signal A: 

P P c(n> 1) - l-c- A = R T . (16) 
This leads to an average number of created photoelectrons for each event with an initial laser signal 

B = A = zMlz*>. (17) 



3.2.1 Determination of the expected number of afterpulses 

For this analysis method all subsequent signals in the time range from 0.1 [is to 16 |is after an initial 
signal are taken into account. After a PMT pulse the baseline might fluctuate causing a fake signal. 
To exclude this effect as well as late-pulses, which appear up to 80 ns following the initial pulse [7 , 
a dead time of 0.1 (J.s after each initial signal is applied. The dark noise distribution is assumed 
to be flat. However, the observed histogram of supposed dark noise signals Bi in figure [3a| shows 
a nonconstant contribution. This is caused by the inefficiency of the setup: signals with a small 
amplitude, which do not pass any discriminator, also produce afterpulses. These afterpulses are 
correlated in time with the strobe signal and contaminate the measured dark noise distribution. 
In order to bypass this effect the average number of dark noise signals B was calculated in the flat 
region from 12 (is to 20 |xs and subtracted from the binned signal histogram Si. For the expected 
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Figure 3: (a) Average temporal distribution of expected dark noise signals Bi for all PMTs. The vertical line 
indicates the flat dark noise contribution B. (b) Average temporal distribution of expected afterpulses Si for all 
PMTs. 



number of afterpulses per SPE event we yield 

£ < = t^f. < 18 > 

with the total number of initial signals N and photoelectrons per initial signal n. The resulting 
average temporal distribution of expected afterpulses of all PMTs is shown in figure |3b| 

Using £i, the expected number afterpulses in the time range up to 16 |J.s after an initial SPE 
signal is given by 

16 \xs 

£ A = E 5 - ( 19 ) 

The average number of expected afterpulses for all PMTs is 

£ A = 2.68 • 1CT 2 , (20) 

with a standard deviation of 

a = 1.46 ■ 1CT 2 . (21) 

To determine the statistical uncertainty the bin entries of Si and Bi are taken to be Poisson dis- 
tributed with an uncertainty of y/~N. The sta tisti cal unce rtainty of n is assumed to be negligible. 
Propagating the uncertainties using equation 18 and 19 yields the statistical uncertainty of £ A . 



As the statistical uncertainty is dominated by the number of measured afterpulses, the relative 
statistical uncertainty (cr s tat/£ A ) is calculated for each individual PMT to allow for a better com- 
parison among the different PMTs. The relative statistical uncertainty averaged over the sum of 

473 PMTs is 

^tlt=(^)=5.1%. (22) 

The only systematic effect is caused by the efficiency e of the setup. It is defined as the ratio of the 
number of measured signals 7V m and the true number of signals iV t . Assuming the same afterpulse 
probability for detected and missed signals, the detection efficiency is given by: 

N N 

N t JV tl ap 
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Figure 4: (a) Average probability distribution Pdj = (IIt=i ^d,fc) ' Pd.i of measuring the first subsequent dark 
noise signal in bin i for all PMTs starting at to = 12 us after the laser trigger. The line is the assumed dark noise 
contribution defined in equation |27| (b) Average temporal afterpulse probability distribution P aPt i of all PMTs per 
bin width At = 160 ns. 



-Wm.ap and Ar ti a P are the number of afterpulses created by the detected signals and the proper 
number of si gnals, respectively. The efficiency can be calculated by usi ng th e the maximum max(£i) 

serves as an estimate for 



of figure 3b as an estimate for AT mjap , whereas max(B i — B) of figure 



the number of missed events. We thus yield 



3 a 



e= m&m) ^ = #- = 0.9991. 

max(£j) + max(i?i — B) 5.505 



(24) 



Compared to the statistical uncertainty this effect can be neglected. 



3.2.2 Determination of the afterpulse probability 

For the second analysis the distribution Si is created containing only the first subsequent signal. 
Dividing the number of signals Si by the total number of initial signals Nq yields the probability 
Ps,i of measuring the first subsequent signal in time bin i. Pg,i is the product of the probability 
that no signal occurred in the time period before bin i, and the probability Pj of measuring at 
least one signal in bin i: 

Ps,i = ^ = l[(l-Pk)-Pi- (25) 

k=l 

Distinguishing between afterpulses and dark noise signals the probability Pg,i becomes 

i-l i-l 

k = l k=l 
i-l i-l 

= (Y[P^,k)-(l[PZ,k)-[Ci--Pd,i)P ap ,i + Pd,i\- (26) 

k = l k=l 

Here, P aPj i and Pa,, are the probabilities of measuring at least one first order afterpulse signal and 
at least one dark noise signal in time bin i. P aPi fe and P^k are the probabilities of measuring no 
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afterpulse signal and no dark noise signal in each time bin k. 

The dark noise probability is assumed to be Poisson distributed with a mean noise rate R: 

n^^e-***- 1 ^*, (27) 
fc=i 

P iti = 1 - e- R At . (28) 

As described in section |3.2.1| the measured dark noise distribution also contains afterpulses due 
to an efficiency of e < 1. For the determination of the dark noise rate R the time range between 
12 us and 20 us after the laser trigger is considered. Counting the number of dark noise cycles D 
containing at least one signal between 12 us and 20 us and dividing it by the total number of dark 
noise cycles Dq leads to the probability of measuring at least one dark noise signal within 8 us: 

D =1 _ e -R-8»s (29) 
Do 

Using this equation the rate R is determined for each PMT individually. 

The probability P ap .i is calculated in an iterative procedure: Starting with the first bin, where the 
probability of measuring no previous signal is equal to one, the probability P ap ,i is given by 

Pap.l = (P S ,l + e- R - w ™ - 1) • e fl ' 16ns . (30) 

Using this result the probability of the following bin Pap, 2 can be calculated. Repeating this 
procedure up to t = 16 us the probability P aPi i is calculated for each time bin. 
To yield an afterpulse probability per SPE event, the bin- wise probabilities are corrected by the 
average number of initial photoelectrons n using equation [4] 

^•^-(l-Pap^n)) 17 ". (31) 

The average temporal afterpulse probability distribution for all PMTs is shown in figure |4b| 
Since the quantities P aPi i are bin-wise probabilities, the sum over several bins does not represent 
any physical quantity anymore. However, the afterpulse probability "P can be calculated by taking 
the advantage of 

16 \\s 

11(1-^) (32) 

4=0 

being the probability of measuring no afterpulse event in the time range from to 16 us. For the 
afterpulse probability of measuring at least one afterpulse in this time interval we can conclude 
that 

16 \xs 

V = l- iJ(l-Pi). (33) 

Figure [5b] shows the total afterpulse probability "P of all PMTs in a histogram. For the average 
afterpulse probability we yield 

P = 2.27-10~ 2 (34) 

with a standard deviation of 

cr = 1.22 • 10~ 2 . (35) 

Since "P is determined in an iterative procedure, the statistical uncertainty is calculated for each 
time bin individually, while the errors of the preceding bins are propagated to the next bin. The 
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Figure 5: (a) Distribution of the expected number of afterpulses S A (equation |l9| for all PMTs. (b) Distribution 
of the afterpulse probability T (equation |33| for all PMTs. 



statistical uncertainty of the entries of the initial signal and dark noise signal distribution are as- 
sumed to be Poisson distributed with an uncertainty of \/N. The statistical uncertainty of the 
average number of photoelectrons n per initial signal is again assumed to be negligible. Propagat- 



ing these uncertainties via equation 26 and computing the squared sum over all a s t a t,i yields the 
statistical uncertainty of "P. As the statistical uncertainty is dominated by the number of after- 
pulses, the relative statistical uncertainty (a s t a t/'P) is calculated for each PMT, which simplifies 
the comparison of different PMTs' results. The relative statistical uncertainty averaged over the 
sum of all PMTs is 

^Lt=(^r)=5.4%. (36) 

Although each uncertainty cr stat ^ is propagated from bin to bin , th e achieved uncertainty is only 
slightly lager compared to the relative error of S A (cf. equation 22). 



The only systematic effect is caused by the efficiency of the setup, as e < 1. Due to high complexity 
of the calculation, impact of this effect on our result was not computed. Nevertheless, we can 
assume that the effect's size is of the same order as in equation |24j and therefore neglect the 
systematic uncertainty. 



3.2.3 Connection between both analyses' results 

The values of S A as well as "P are smaller than 0.1 for all PMTs, as shown in figures 5a and 5b This 
means not only the expected number of afterpulses per SPE is less than 0.1, but also probability 
for a single PE to obtain at least one afterpulse is smaller than 10%. Therefore the selection 
specification of the Double Chooz experiment were well fulfilled, which requested an afterpulse 
probability of less than 10 % in the time range of to 16 us. 



The correlation between £ A and T is plotted in figure [6a| In figure 6b we can see that the 
temporal distribution (gained via "Pi together with equation [4] and equation [3]) approaches zero 
for t > 12 [xs, i.e. no further first order afterpulses occur after 12 us. Thus, the calculated probability 
V includes all created first order afterpulses. On the other hand, the temporal distribution Si does 
not approach zero for times larger than 16 us. Therefore the expected number of afterpulses £ A is 
calculated by summing Si from to 16 us, including only a fraction / of higher order afterpulses. 
This observation shows, that the fraction / has to be considered in our calculations. 
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Figure 6: (a) E A vs. "P for all PMTs. (b) Temporal distributions of the expected number of afterpulses S A and the 
afterpulse probability fj,i, averaged over the sum of 473 PMTs. 



According to equations [3] and 13 
can be calculated using £ A and "P: 



the average number of photoelectrons I contained by an afterpulse 
£ A + ln(l-7>) 



(37) 



£ A • ln(l -V) + (ln(l - V)) 2 ■ (1 - /) ' 

The fraction / is calculated using the temporal distribution j of expected higher order afterpulses 
(cf. equation 151. To determine E^i via equation 14 we use the average temporal distributions of 
Si and /Ltj (cf . figure 6b I . This leads to an average fraction of 

7=82.6%. (38) 

Compared to the uncertainties of the other values used to calculate £, the uncertainty of / can be 
neglected. 

Using /, £ A and "P, 1 is calculated for each PMT individually. For the error weighted mean of all 
PMTs we receive 

2= (4.95 ±0.13) PE. (39) 



4 Timing and charge measurements for dedicated PMTs 

In order to determine the temporal occurrence of afterpulses correlated with their charge, an 
additional and more detailed measurement was realized. Due to complexity these tests were 
performed only for a few PMTs. To gain information about the dependency of the afterpulse 
properties on the incident light levels, different intensities were used. 



4.1 Experimental Setup 

The tested PMTs were arranged in a light-tight Faraday room as described above. All electronics 
components were installed outside the room, the setup is shown in figure [7j A custom-built 
LED board with a wavelength of 380 nm served as light source. To allow for a simultaneous 
measurement of timing and charge of the afterpulses, both a time-to-digital converter (TDC, model 
V775/ CAEN) and a charge-to-digital converter (QDC, model V792/ CAEN) were used. Another 
VME module producing rectangular pulses with a distinct frequency, acted as pulse generator. 
It provided several output channels and the pulses of each output could be adjusted individually 
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Figure 7: Scheme of the experimental setup for detailed measurements. 



in width and in their delay with respect to the other signals. In combination with discriminator 
modules and a logical circuit composed of AND / OR and gate delay modules, pulses occurring after 
an initial PMT pulse could be measured. 

The TDC used for this measurement was a single-hit TDC with a resolution of 0.3 ns per bin. 
It records the first appearing afterpulse in a time window of approximately 1200 ns. Therefore 
several measurements at delay steps of 1 us were carried out. Due to baseline distortion effects 
caused by the initial PMT pulses, the search for afterpulses was started 375 ns after an initial pulse. 
Assuming that the dark noise is independent in time, the measured time and charge spectra of 
the afterpulses are corrected by taking an additional dark rate measurement with a time delay of 
~ 70 lis with respect to the initial pulse. 



4.2 Analysis and Results 

The results presented in this section characterize the afterpulse properties of one particular PMT 
in the range A = [375 ns; 12 us] after the main PMT pulses, as the afterpulse behaviour turned out 
to be similar for all tested PMTs. 

Due to usage of a single-hit TDC several individual measurements covering different regions after 
a main pulse have been taken. Merging the datasets results in steps at each edge to a subsequent 
measurement range. This happens since events in earlier bins suppress possible afterpulse events 
in any of the subsequent bins of one independent measurement. These kinks can be removed by 
the following bin-by-bin correction of the spectra. The correction assumes that the probabilities to 
obtain afterpulses in certain bins do not depend on each other and the maximum expected num- 
ber of afterpulses per bin is smaller than 1. This is implied by the appropriate assumption that 
multiple afterpulse events produced by a single photoelectron are independent events and that 
the contribution of higher order afterpulses within a measurement interval of 1 us is negligible. 
Coincidental afterpulses in the same bin can be excluded due to the low afterpulse probability in 
a bin width of 0.3 ns. 

The first bin of a dataset does not need to be corrected. The corrected number of entries for the 
second bin is given by 

ocorr ^2 

" 2 - i-sm 

where No is the number of initial pulses. The correction for higher bins accounts for the capability 
of one or more events in any prior bin to suppress afterpulse events in the bin of interest. Hence, 



(40) 
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the number of events of the i-th bin Sf on ' is 



sr r = $- n^-^ "/^) • (4i) 



\k=l 



Figure [8a]shows both the measured temporal afterpulse spectrum in black and the corrected tempo- 
ral spectrum in red. The additional dark rate measurement has as well been corrected via equation 
and was used to calculate the mean number of dark rate signals per bin B. The spectrum Sf°" 
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was dark noise corrected in each bin by subtraction of B. 



The histogram with all time differences between initial pulse and afterpulses is plotted in figure 8b 
representing the afterpulse timing spectrum for the measurement with n = 1.1 PE. Distinctive 
peaks can be found at around 2, 6.5 and 8 [is. Comparing the afterpulse timing spectra, the trigger 



board (figure 6b I and the TDC measurement are in agreement for times greater than 4 us, for 
smaller times the TDC measurement shows an additional exponential contribution. To investigate 
this observation further, a control measurement using the setup shown in figure [7] with a differ- 
ent type of PMT was performed. As expected a differing peak structure but also an identical 
exponentially decreasing part was observed in the timing spectrum. Furthermore, measurements 
using different LED intensities could demonstrate a non-linear dependence between the exponential 
component and the LED intensity. This led to the conclusion that the exponential contribution is 
induced by afterglowing effects of the LED. 



The afterpulse charge spectrum (figure 9b and 9b I is dominated by small charges of a few photo- 



electrons. However, a contribution of higher charged afterpulses up to 40 PE also becomes evident. 
The exponential tail in the timing due to LED afterglow is associated with charges £ < 2PE, as 
shown in figure [9c] Also, afterpulses with a timing of 5-11 us contain a considerable amount of low 
charges. Regions of afterpulses with higher charges do not fully correspond to the peaks seen in 
the time spectrum in figure [85] In figure [9a| and [9c] regions of afterpulses transferring large charges 
can be identified at 500 ns, from 1.5 to 2 us and around 8.5 |xs. 

The effect of suppressing possible subsequent pulses due to usage of a single-hit TDC can not be 



corrected in case of the charge spectra (cf. equation 41 for the time spectra). Hence, a slight 
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deformation of each charge spectrum is expected. 

Measurements with higher light levels confirmed the assumption that highly charged afterpulses 
arise due to ion feedback. If more light is focused on the photocathode, and thus a larger number 
of initial photoelectrons n is released, the number of afterpulses rise, but the spectral shape of 
timing and charge do not change. 

Furthermore, additional afterpulse measurements with three different voltages were performed 
for the same PMT for n = 14.8 PE. The conventional high voltage HVi leads to a gain of 0.8 pC 
per photoelectron, HV2 provides a factor two in gain (1.6pC per photoelectron) , and HV0.5 results 
in half the gain of HVi. No significant shift of the peaks structure could be observed. The change 
in high voltage showed only a slight effect on the number of afterpulses as seen in figure |9d| 

5 Discussion of the shape of the temporal afterpulse distribution 

As explained in section [2] afterpulses are caused by ionization processes between the photocathode 
and the first dynode. The time distribution of the afterpulses is dominated by the rather slow 
movement of the heavy ions with respect to the transit time of the photoelectrons. Due to the 
nonuniform electric field of a hemispherical photomultiplier the travel time between the ionization 
location and the photocathode is not sensitive to the ionization location [S]. The travel time can 
be calculated to 

< 42 » 

Here, m and I are the mass and charge of the ion and Vo is the voltage between photocathode and 
first dynode. 



Both temporal afterpulse distributions £ t and \ii show the same peak structure (cf. figure 6b I 
with at least five peaks with mean times approximately at 0.52 us, 2.1 ujs, 4.9 us, 6.3 \is and 7.9 |xs. 
This leads to the assumption that the afterpulses are mainly caused by five different types of ions. 
According to calculations from Hamamatsu the peak at 2.1 us is caused by Methane (CH4) and the 
peak structure between 6-8 us is caused by Caesium (Cs). Other possible candidates are Hydrogen 
(H) and Helium (He) as well as the cathode materials Potassium (K) and Antimony (Sb). However, 
it's not possible to clearly assign these elements to the observed peaks. 

6 Conclusion 

We presented afterpulse timing measurements for 473 photomultipliers which served as one of the 
qualification tests for the Double Chooz experiment. The selection specification of the apfterpulse 
probability has been limited to 10% for each PMT. It was shown that all PMTs fulfilled this 
requirement (figure [5]). 

The temporal distribution of the afterpulse occurence was determined for the number of expected 
afterpulses. In a second analysis using the same data set the probability of measuring at least 
one afterpulse and the timing spectrum for the expected number of first order afterpulses was 



computed. Both methods' results showed as expected a strong correlation (cf. figure 6a) and were 
used to derive the average charge transferred by an afterpulse I. 

More detailed timing and charge measurements of afterpulses were performed for a few PMTs 
using a different setup. The results of one PMT was presented, which represents well the general 
behaviour (figure |9| . Comparing the afterpulse timing spectra, both measurements show a good 
agreement for times greater than 4 us. The additional exponential contribution in the second 
measurement has been identified to be caused by afterglowing effects of the LED. 
The two different analysis methods used for the first measurement allowed us to determine the 
average number of photoelectrons contained by afterpulses of £ = (4.95 ± 0.13) PE. As the charge 
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Figure 9: (a) Scatter plot (t,£), temporal occurrence t vs. the charge I of the afterpulses for n = 17.4 PE. (b) 
Charge spectrum of the afterpulses for n = 1.1 PE. (c) Afterpulse timing spectrum with charge cuts for n = 1.1 PE; 
black: complete spectrum, green: < £ < 2PE, blue: 2 < I < 10 PE, red: 10 < £ < 25 PE. (d) Afterpulse timing 
spectrum for n = 14.8 PE measured with different high voltages: blue: HV0.5 = 1130 V, black: HVi = 1229 V, red: 
HV 2 = 1340 V 
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measurements contain also LED afterglowing events, a direct comparison of the average number of 
PE carried by afterpulses between the first and the second measurement was not possible. However, 
the charge measurements indicate that afterpulses can carry charges in the range of 1 and 40 PE 
(figure 9a and 9b I . 

The measurements performed with different high voltages does not show a significant effect on 
the afterpulse time distribution shape. Only a small linear correlation between the applied high 
voltage and the afterpulse probability can be observed. 
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